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^ ■ Abstract 

£^ ' A general formula for the linearized Poincare map of a billiard with a potential is 

\q . derived. The stability of periodic orbits is given by the trace of a product of matrices 

describing the piecewise free motion between reflections and the contributions from the 
reflections alone. For the case without potential this reduces to well known formulas. Four 
billiards with potentials for which the free motion is integrable are treated as examples: 
The linear gravitational potential, the constant magnetic field, the harmonic potential, 
and a billiard in a rotating frame of reference, imitating the restricted three body problem. 
The linear stability of periodic orbits with period one and two is analyzed with the help 
of stability diagrams, showing the essential parameter dependence of the residue of the 
periodic orbits for these examples. 

Q\ '. PACS: 03.20.+i, 05.45. +b 

la: 

■ 1 Introduction 

In the theory of dynamical systems billiards were introduced in the twenties of this century by 

■ Birkhoff. His idea was to have a simple class of models that show the complicated behavior 
of nonintegrable smooth Hamiltonian systems without the need to integrate a differential 
equation [1]: "But in this problem the formal side, usually so formidable in dynamics, almost 
completely disappears, and only the interesting qualitative questions need to be considered" . 
Birkhoff studied billiards inside smooth convex domains, for which it is now known [2] that 
they can not be ergodic. In the seventies Sinai proved the ergodicity of the billiard on a 
torus with a disk removed [3], and subsequently Bunimovich, Wojtkowski and others [4, 5] 
designed larger classes of ergodic billiards. Starting with some early papers in the eighties 
[6, 7] billiards became quite popular in physics. For investigations of both, the classical 
dynamics and the quantum mechanical properties, periodic orbits and their properties are of 
major interest, which are investigated using the Poincare map. 

For billiards the natural Poincare map is from bounce to bounce: The boundary <90 of 
the billiard domain C M 2 is parametrized by the arc length s. d£l is assumed to be C 2 , 
such that the curvature is well defined. Reflections at points in the boundary which are less 
smooth need a special treatment. Starting at some point r(si) £ <90 all possible orbits leaving 
<9f2 can be parametrized by an angle (3 from —tt/2 to tt/2, where —tt/2, 0, tt/2 specify the 
directions — t(si), n(si), t(si) where t(si) is the tangent vector to <90 at r(si) and n(si) is 
the normal vector pointing into O. We have d s r(s) = t(s) and d s t(s) = re(s)n(s) where k(s) 
is the signed curvature at s, negative k indicating a defocusing boundary with respect to O, 
or, equivalently, the normal n. If instead of (3 the component p of the velocity parallel to t is 
used, p = sin/3, then (s,p) are canonical coordinates on the Poincare surface of section. The 
trajectory starting with £i = (s\,pi) encounters the boundary after time r at S2 and has a 
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parallel component P2 (with respect to t(s2)) before and after the reflection. The Poincare 
map is denoted by P such that P£i = £2- 

A periodic orbit of period n is now described by a list of points £j = (si,pi), i = 1, . . . , n 
such that P^i = for i = 1, . . . , n — 1 and -P£ n = £i, or equivalently P n £,\ = £i. Finding 
periodic orbits is a hard problem, even for billiards. In this paper we take the point of view 
that the periodic orbits are given to us, and we want to find a formula for their stability, i.e. 
for the linearized map DP of P. 

In the following we will denote a quantity to be evaluated at £j or by a lower index i, 
e.g. r(sj) = r-j, t(s«) = tj or DPj^i for the linearized map of P from £j to Then we 

obtain for Z)P n 

5^,1 = £^n,n-l • • • DP 32 DP 21 (1) 

by the chain rule. For a smooth Hamiltonian system the linearized map depends on the solu- 
tion of the variational differential equation, which can usually only be obtained by numerical 
integration. For an ordinary billiard without potential, rather than the infinite amount of 
information contained in this solution, we just need three numbers per reflection, besides the 
initial and final point: the time r (or length) between reflections, and the curvatures K{ and 
Ki+i of d£l at the initial and final points and r^_|_i . The expression of DP in terms of these 
data is well known, see e.g. [8, 9, 10]. The main point in stability formulas of this type is 
that they do not express everything in terms of the initial point £, which is of course possible. 
Instead, some quantities from the free motion, e.g. the time of flight r, directly enter the 
formula. The time r is usually given by the solution of a complicated equation (only for 
special boundaries can it be obtained from the roots of a polynomial of degree up to four). 
Therefore it makes sense (also numerically) to use r in stability formulas, because it must 
have been already calculated to obtain P. 

If we add a simple potential to the billiard, the problem of the calculation of DP becomes 
a little harder. To keep the original idea of Birkhoff it only makes sense to add integrable 
potentials for which the solutions of the free motion in the potential can be written down 
explicitly. The general formula we obtain does not depend on the integrability, however 
in the examples we only treat integrable potentials. Typical examples are linear potentials 
[11, 12, 13, 14, 15, 9, 16] with free motion along parabolas, quadratic potentials [17, 9] with free 
harmonic motion, magnetic fields [18, 19] with free motion on circles and the fictitious forces of 
a rotating frame of reference [20]. The reason for the study of billiards with potential, besides 
of course the situation where they model a given physical situation, e.g. the restricted three 
body problem in [20] , is that the potential can introduce two typical features of Hamiltonian 
systems which are absent in ordinary billiards. The most important is the energy dependence: 
As it is nicely illustrated in [15], one can observe the typical evolution from an integrable low 
energy limit, via an intermediate regime with the most complicated motion to an integrable 
high energy limit. This scenario is typical for bounded Hamiltonian systems, and billiards 
with potential are simple model systems of this kind. The second reason is that a magnetic 
field (or a rotating system) breaks the time reversal symmetry which is important because it 
influences the bifurcations of periodic orbits. 

The paper is organized as follows. In sec. 2 we take a global point of view and discuss 
the topology of the surfaces of section induced by the potential. The linearized Poincare map 
is obtained in sec. 3, and the results are used to calculate the stability of periodic orbits in 
sec. 4. In the remaining sections we calculate the stability of short periodic orbits for billiards 
with linear (sec. 5), quadratic or magnetic potential (sec. 6), and a rotating system, which 
models the restricted three body problem (sec. 7). 
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2 Billiards with potential 



Let the Hamiltonian of the free motion be given by 

^(q,p) = ^(p-^(q)) 2 + ^(q). (2) 

The velocity is given by v = p — A. In principle the calculation could be done with an 
additional positive definite "metric" in the kinetic energy, but in the application we have in 
mind we only need the above form. Let the configuration space of the free motion be M 2 . 
The energy surface of the billiard with potential is just the ordinary energy surface but with 
q restricted to the interior of the billiard Q, 

4 = {(q,p)|%p)=/i,qefi}. (3) 

Let the solution of Hamilton's equations X = JVA be given by the flow <\> 

<l> t :X 1 ~X(X 1 ;t), (4) 

where we introduced X = (q, p) and the symplectic matrix J. The section condition in our 
system is naturally given by the billiard boundary 3f2 C M 2 or explicitly by q(t) = r(s), such 
that the surface of section becomes 

£ ft = {(q,p)|fT(q,p) = h,q€dSl}. (5) 

On the surface of section T,^ we introduce coordinates £ = (s,p), where s is the arc length of 
dfl, and p = (t, v) is the component of the velocity v = p — ^4(q(s)) parallel to the tangent 
vector t(s). 1 As we will see in the next section this is a canonical coordinate system on 
Y<h and therefore P is area preserving in these coordinates. It is convenient to extend this 
coordinate system on £/,, to the neighborhood of S/j in phase space. For this we take as local 
coordinates the energy h in the direction of X7H, and the time \x in the direction of the flow 
JVH. The coordinate tuple (s,n,p, h) will be denoted by S. 

The Poincare map P from bounce to bounce in a billiard with potential can now be 
described by a succession of four maps: 

the map from the surface of section to phase space, Hi = (s\, 0,p±, h) 
the flow map of free motion in the potential 
the map back to local coordinates on the surface of section 
the reflection, E 2 = (s 2 ,0,p2,h) 

(6) 

The time t\ must be chosen in such a way that X2 £ S/j, where the tilde denotes the moment 
before reflection, see below. As we will see, the description as a four dimensional mapping is 
convenient for the linearization of the map. If we write a£ we really mean aE where S is the 
trivial 4-D extension of £, similarly for er -1 . 
In this notation the Poincare map reads 

P : H h — > T, h 

i 1 ^ P(£) = P a~V T ^; (7) 
with r determined by 4> T o(, £ S^. 

Note that if there is no reflection (i.e. if we study a Poincare section of the free Hamiltonian) 
er -1 from the last iteration and a from the next cancel, such that the iterated map is given 
by 

P n ^ = (j ^ (p Tn ■ ■ ■ tfPafr = <j-V {T1+ "' +Tn) ^i (8) 
1 Note that the p stands for parallel component of the velocity, not of the momentum p 



a : Hi ^ X\ 

(f) Tl : Xi h+ X 2 

cr _1 : X 2 i-> H 2 

p : ~ 2 >-> H 2 
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as must be the case. 

The tilde on X2 and H2 indicates that we are on the surface of section but still in the 
instant before the reflection, such that the normal component n,2 = (112, V2) of V2 is smaller 
than 0. Here we take the point of view described in [21] that the surface of section consists of 
both parts with either sign of n<i- In the case of an ordinary billiard with a smooth boundary 
Y<h becomes a 2-torus. This 2-torus is divided into two cylinders by the lines where the 
flow is tangent to the surface of section. They are given by the condition that the velocity 
v is parallel (or anti parallel) to the tangent t, or in local t-n-coordinates by n = 0. In 
the case of a convex billiard these lines are two periodic orbits given by the clockwise and 
counterclockwise sliding motion along the boundary. One cylinder carries initial conditions 
with the velocity pointing inward (i.e. after the reflection, n > 0) while on the other one 
the velocity is pointing outward (i.e. before the reflection, n < 0). The flow 4> T maps the 
first cylinder to the second, and the reflection p maps the second back to the first. In local 
t-n-coordinates p is trivial: it just changes h to —n. This is the reason why we first use a 1 
and then p. Actually the points on S/j with outward pointing velocities can be considered the 
same points as the ones with the reflected inward pointing velocity, such that the reflection 
map gives an identification of the boundary of the filled torus ft x S 1 , which is the energy 
surface, and turns Sh into a manifold without boundary. In the present case we obtain S 3 by 
this construction. 

If we add a potential the topological situation changes: the accessible region of motion in 
configuration space Qh can now be a subset of the region Q bounded by the billiard walls. 
As usual Qh is given by the condition U(q) < h, such that for the billiard with potential we 
obtain 

Qh = {ci\U(c l ) <M^}. (9) 

The additional boundaries dQh\dd are the ovals of zero velocity [1], where the kinetic energy 
vanishes due to the additional potential. Moreover there might be a subset of the energy 
surface, such that orbits in this subset never meet the boundary of Q. In the examples we are 
only dealing with integrable free Hamiltonians such that this part of the energy surface will 
be foliated by tori. Although the Poincare section will be ^-incomplete [21], i.e. not every 
orbit in the energy surface will intersect S^, the missing orbits are well under control and 
we ignore them in the following discussion. A very important point in adding a potential is 
that Shi the dynamics, and also T,h in general depend on the energy: instead of the entire 
<9$7 the part of the boundary in the accessible region of configuration space can now become 
a collection of intervals / C dQ if the energy is not too high. On every point of / we have 
a certain amount of kinetic energy available such that the norm of the velocity is fixed and 
all possible directions can be parametrized by an angle, such that we have to attach circles 
S 1 to every point of I. However, at the boundaries of / where the ovals of zero velocity 
intersect the billiard boundary dQ, this circle has to be contracted to a point since the only 
possible velocity is 0. Therefore the surface of section becomes a collection of spheres, one for 
each interval in /. On each sphere there is a line where the velocity is tangent to d£l, and it 
divides the sphere into two disks. If we restrict our attention to the moment after reflection, 
the Poincare map will be defined on this collection of disks. 

3 The linearized map 

For a general discussion of a Poincare map from a (curved) surface in the energy surface to 
itself see e.g. [22, 23, 24]. As we will see, a symplectic coordinate system in the neighborhood 
of Eft is given by the arc length s, the parallel component p of the velocity v (not of the 
momentum p), the value of the Hamiltonian h, and the flow time p starting at the point 
specified by s,p, h. Denote the components of this coordinate system by S = (s, p,p, h). The 
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Figure 1: The orbit q(i) encounters the boundary dil given by the arc length parametrization r(s) 
at r with velocity v before the reflection. The angle between the tangent t at r and a fixed horizontal 
reference direction is a. The angle between the reflected velocity v and the normal vector n is (3. 
In the t-n-coordinate system the velocity v is given by (p, n), such that p — |v| cos/3. Usually all 
quantities have an index to denote the point at which they are to be evaluated. 



transformation to the original canonical variables X = (q, p) is given by 

X = aZ = <ffF(8,p,h), (10) 

where F gives the initial point on in the coordinates X. The position is given by q = r(s) 
and the velocity is v = D(a(s))u where u = (p, n) is the velocity in the local t-n-coordinate 
system and D = (t,n) the rotation matrix, a is the angle between t and the horizontal 
axis. To fix an orientation of the t-n-system we let n always point into Q and then rotate 
by S = D(7r/2) to obtain — t = Sn. By this convention the curve must be traversed 
counterclockwise for the billiard inside the curve and clockwise for the billiard outside the 
curve. This automatically fixes the sign of the curvature k to be positive for the billiard 
inside a convex domain f2, and negative for the corresponding outer billiard. 

The normal component n is implicitly defined via the shifted energy function on given 

by 

E(s,p, n, h) = \{P 2 + n 2 ) + U(r(s)) - h, (11) 

where the positive square root gives n and the negative n. Now F is given by 

F:(s,p,h)^X = (q,p)* = (r(s), D(a(a))u + A)'. (12) 

We only want to know the transformation (10) close to S^, i.e. for small [i such that we can 
approximate ^ ~ 1 + fiJVH, where the Hamiltonian vector field X = JX7H at r(s) in the 
coordinates X is needed, 

^(^""'((W-v)' (13) 

The complete transformation, which gives a for ji = 0, is given by 

X = F(s,p,h)+»X\ F{s ^ h) =( l + Hy (14) 
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with p and p given by (12,13) and n implicitly determined by (11). Now we want to find the 
linearization of this map at jjl = and show that it is a symplectic mapping. By this we find 
that the variables £ = (s,p) are canonical on S^. Moreover we need exactly this symplectic 
mapping to determine the linearized Poincare map. Using (11) we obtain the derivatives of 
n. With the abbreviation 7 = (d q U,t) we find 

dn 8E ,dE 

~5T = ~Pi -U/ n - ( 15 ) 



d(s,p, h) d(s,p, h) I dn 
The Jacobian of the transformation (14) at u = is the 4 x 4-matrix 



ax 
dz 



= (d s F,X,d p F,d h F) (16) 

tv \ _ / A; 
d s p p -Sv/n n/n / ' I Q D; 



(17) 



S = D(tt/2), and with (12) we find 



d sP = (d s D)u + m s u + {d q A,d s r) (is) 

= kSv -n<y/n + (d q A,t). (19) 

The symplectic condition, L*JL = J, states that A^D; = 1 and that AjCi be symmetric. 
The first holds because of t*n = 0, — v*Sv = and v*n = — t*Sv = n, such that D; = Af u 
and det A^ = n. The second condition is the identity 

<v,a iP ) = (t, P ) (20) 

-(v,n) 7 /n + (v,^t) = (t,a,^*v) - (t,d q U) (21) 
- 7 + <v,a,At) = {t,d q A t w)- 1 . (22) 

The inverse of L can be simply computed because L is symplectic: 



-1$ vO ^={%^- 

To linearize the Poincare map we use L|^ = Li to map vectors at £1 from TS^ to Euclidean 
coordinates, which are then mapped by the Jacobian of the flow M = deft /dX. At the new 
point X2 we return to the section coordinates with L _1 and keep in mind that we must still 
perform a reflection. In the local coordinate system this just means to put — 112 instead of fi2 
in L^" 1 . To indicate this we write Lj = L(si,pi,rii) = h(si,pi,hi) = L(si,Pi, — rij) such that 

a — 1 

-^L^MxLl (24) 

To obtain a 2 x 2 matrix for DP we just set dr and dh equal to zero, i.e. we consider the 
restriction of the map to the tangent space TT,^ in local coordinates on E^. Deleting the 
appropriate rows and columns from L and L _1 we find 

where the left matrix is a 2 x 4-matrix and the right is a 4 x 2-matrix. For the trivial 
Hamiltonian of the ordinary billiard the only term in d s p is related to the curvature, see 
(19). In the presence of a potential U there is an additional contribution in the form of 
7 = (d q U,t), and a similar term (d q A, t) for the vector potential. 
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The Jacobian for the free motion in ordinary billiards is 

M=(J7), P6) 

and the main simplifications in this case come from the fact that V2 = Vi, and with (19) we 
can rewrite (25) as 

DP21 = ( -K 2 n 2 1 ) ( ^O 7 " 12 t" ) ( ^ ) ( -Svi/m ) ( -lm I 2 } 

-(*:)(-*?)(? =ou !)(oi) « 

which is a well known formula, see e.g. [8, 10]. For cases with potential more explicit formulas 
than (25) tend to be quite complicated. As we will see, it is nevertheless possible in some 
cases to perform the stability analysis up to period two. The obvious use of (25) is in the 
numerical calculation of Lyapunov exponents (see e.g. [10, 25]) and in the numerical search 
for periodic orbits with Newton's method. Note that if different, possibly non-canonical, 
coordinates are used, as the arc length parametrization of d£l may be too complicated, the 
expression for DP21 has to be sandwiched by the Jacobian of this transformation for Newton's 
method to work. Another application is the analytical estimation of the Lyapunov exponent 
as performed in [25]. Let us remark that the linearized map in the full phase space (24), 
besides being conceptually much clearer than (25), can be useful to implement a continuation 
method to find orbit cylinders in phase space, i.e. the collection of periodic orbits for different 
h, because it contains information about the change of the orbit if the energy changes. 



4 Linear stability of periodic orbits 

The standard application of the linearized map is to calculate the stability of periodic orbits. 
There is one simplification in calculating the stability of periodic orbits as compared to the 
calculation of the general linearized Poincare map, namely that we are only interested in the 
trace of a product of linearized maps, and that we can cyclically permute matrices under the 
trace. The 4x4 monodromy matrix of a period n orbit is given by 

and using (24) we obtain 

trW := trM£) = tr (l^M^ L^M^L^ • • • L^M^) . (30) 
Now we define the linearized reflection matrix 

K, = L-%, (31) 

and with a cyclic permutation we find 

n 

trW=trJjM i K i . (32) 
i=i 

The stability is split up into contributions from the reflection K and from the free motion 
M. Moreover, both contributions depend only on a single point, while for we needed 
information about the initial and final point. A similar formula with 2x2 matrices for the 
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ordinary billiard is well known, see e.g. [26, 27]. Multiplying DP32DP21 as given by (28) we 
find 

? ) ( 2re 2 /n 2 -1 ) - (33) 

In (32) we have the natural generalization for an arbitrary billiard with potential, which is 
conveniently formulated in four dimensions. 

The reflection matrix defined in (31) can be written as 

k =(rw)- <«> 

It is symplectic, which amounts to WW* = 1 and W*R symmetric. W is a symmetric 
involutive matrix with det W = —1 and W 2 = 1, as must be the case for a reflection: 

I sin 2a — cos 2a J 

The eigenvalue-eigenvector pairs of W are (+l,t) and (— l,n). The reflection line is t and 
the direction of reflection is n, recall that a is the angle between t and the horizontal axis. 
Note that the reflection W acts on both, momenta and positions. In this way the mapping 
"corrects" the time at which the reflection took place: The symplectic coordinate system H 
ensures that we do not have to bother to directly calculate the (linear approximation of the) 
change of time between reflections. 

The lower left block R of K is much more complicated than the other blocks and contains 
all the contributions from re, U and A. From (31) and (23) we find after some manipulations 

re 2 Ad 

R = 2-SDu(SDW u)* + -(&[/, n)nn* - — Bnn 1 + 4(n, d g At)8, (36) 
n n n 

where Wo is the above reflection matrix at a = and with q = (x, y) we introduce B = 
d x A y — d y A x , the magnetic field. Except for the skew symmetric contribution from d q A the 
matrix R is obviously singular. For constant magnetic fields we have (n, d q At) = B/2. In 
the following we will make this assumption, which is valid for the examples treated below. 
Note that SDW u = SWDu = SWv such that we obtain 

2k 2 2n 

R = — Sv(SWv)' + -(0,17, nW - 2B(-±-nn t + S). (37) 
n n n 

In the stability formulas we will often encounter 

2 

trR = - (k( P 2 - n 2 ) + (d„U, n) - 2pB) , R 21 - R 12 = 4rep - 4B. (38) 
n 

The monodromy matrix has eigenvalue 1 twice, with left eigenvector X7H and right 
eigenvector JX7H [24, 28]. In the S coordinate system it has the special form 



Ml 4 ) 



/ a b x \ 

X 1 X X 

c d x 

\ 1 / 



= Lr 1 M„K n ---M 1 K 1 L 1 , (39) 



where x stands for arbitrary coefficients, and the entries a, b, c, d give the linearized 2x2 
Poincare map DP n in local coordinates. From det = 1 we see once more that det DP = 
ad — be = 1, such that £ are area preserving coordinates for P. Right now we are interested 
in the fact that the spectrum of is {1, 1, A, 1/A} such that 

tr (4) =2 + A + _L. (40) 
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By calculating tr„ we therefore obtain the eigenvalues of DP n : 

tr^ 2) := tr DP n = trM^ - 2 = tr^ - 2 = A + -. (41) 

A 

Period two orbits in the ordinary billiard 

Let us now calculate tr^ for a period two orbit of the billiard without potential, where M 
is given by (26). It will turn out that the same M holds for the case of the gravitational 
billiard. 

t4 4) = tr(M 2 K 2 MiKi) (42) 
= 4cos(2q! - 2a 2 ) + (n + r 2 )tr(W 2 Ri + R 2 Wi) + nr 2 trR 2 Ri (43) 

For a period two orbit (without potential) we have p± = p 2 = 0, n\ = n 2 = 1, r = t\ = r 2 
and a\ = a 2 + it, hence we obtain 

t4 4) = 4-4r(«i + k 2 ) + 4kik 2 t 2 , (44) 

which is a well known result, see e.g. [9]. Before passing to the new application to gravitational 
billiards, we show the stability diagram of this period two orbit. It is convenient to introduce 
Greene's residue R [29] and its "complement" R defined by 

R = 1 - R = tr^/4, (45) 

such that in the above case 

1 - R = R = (1 - r/ei)(l - tk 2 ). (46) 

Bifurcations with A = +1, R = are called "direct parabolic" and bifurcations with 
A = — 1, R = "invers parabolic". The direct parabolic bifurcation with A = 1 generically 
means the creation of periodic orbits "out of nothing" in a saddle-center bifurcation. Thus, 
the continuous dependence of orbits on the bifurcation parameter is not guaranteed at this 
point. Index preservation [30] requires the simultaneous creation of one orbit with index 
+1, which must be elliptic close to R = 0, and one orbit with index —1 which is always 
direct hyperbolic. In the presence of symmetry (e.g. time reversal symmetry) a symmetric 
orbit may depend continuously on the bifurcation parameter, and in this case one observes a 
symmetry breaking pitch fork bifurcation instead [31]. 

At an invers parabolic bifurcation with A = —1, R = the periodic orbit depends contin- 
uously on the bifurcation parameter, and this is the most important difference as compared 
to the direct parabolic case. The orbit has index +1, and there is no way to create an orbit 
with the same period n and index —1 because R = 1. Therefore the orbit must survive the 
crossing of the line R = 1, and turns from elliptic via invers parabolic to inverse hyperbolic. 
However, considering the map P iterated 2n times, an inverse hyperbolic orbit of period n 
looks hyperbolic, such that we have index —1 for P 2n , and therefore we typically expect the 
creation of two stable orbits of period 2n with index +1 (or the destruction of two unstable 
ones). This is called period doubling bifurcation. Of course there can be more complicated 
bifurcations at A = ±1, and also additional resonances in the elliptic region, see e.g. [30, 32]. 

In Fig. 2 we show the level lines of R as a function of tk\ and tk 2 . We could scale the 
length by r such that we always deal with a unit length orbit, but sometimes it is more 
natural to think of an orbit as becoming longer instead of the curvature becoming larger. 
Regions with elliptic orbits (0 < R < 1) are shaded. The boundaries are given by R = (i.e. 
by the hyperbola t = 1/ki + 1/k 2 ), the direct parabolic case, and by the two lines R = 
with tki = 1 or tk 2 = 1, the invers parabolic case. The interpretation is as follows, where 
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-2-10 1 2 3 

XKi 

Figure 2: Stability diagram for period two orbits in the billiard without potential. Contour lines 
of the residue R (46) are shown. The elliptic region is shaded. Inside this region the lines with 
R = 1/2,3/4 are drawn. In the unstable regions the integer valued residues are shown. The bold 
dashed line indicates the stable period two orbit in the family of billiards with elliptic boundary. 
Period two orbits of the Sinai billiard are in the lower left quadrant, those of the stadium billiard, or 
the unstable ones in the ellipse, are in the upper right region with R < 0. With an appropriate change 
of the axes labeling this diagram can also be used for period two orbits of the billiard in a constant 
magnetic field (56) and for period one and period two orbits of the billiard with a harmonic potential 
(65,66,67). 

due to symmetry we restrict the discussion to TK2 > tk\\ 

If both boundaries at the reflection points are defocusing, Kj < 0, the orbit is unstable. The 
direct parabolic reflection at two parallel straight lines is at the upper right corner of this 
quadrant. If there is one focusing and one defocusing reflection the orbit is usually inverse 
hyperbolic; however, even for large defocusing curvature, elliptic motion is always possible. 
If we lower the value of TK2 starting in the inverse hyperbolic region we will encounter an 
inverse parabolic bifurcation at TK2 = 1. Then there is a small elliptic region which we leave 
via a direct parabolic bifurcation. There is no further bifurcation if the boundary becomes 
defocusing. 

The region with two focusing reflections, «j > 0, is the most complicated. If the orbit is 
sufficiently long, then it is direct hyperbolic, R < 0. This is the longer period two orbit that is 
always present in a smooth convex billiard. If we make the orbit shorter we encounter a direct 
parabolic bifurcation if tk\ is sufficiently close to 1. Examples of orbits in the intermediate 
stable region are the shorter period two orbit of smooth convex billiards (if they are stable). 
At tk\ = 1 these orbits become inverse hyperbolic if the curvatures are different, even though 
both reflections are focusing. However if the length is decreased further the orbit becomes 
stable again at TK2 = 1. 

The period two orbit of the circle billiard is at tk\ = TK2 = 2, and for an ellipse with semi 
major axes a and b, a > b the unstable orbit along the longer axis is at TKi = 2a 2 /b 2 while 
the stable one is at TKi = 2b 2 /a 2 , indicated by the bold dashed line. For a = by/2 it is invers 
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parabolic, such that this situation is particularly susceptible to perturbations. Note that the 
passage of the stable periodic orbit of an integrable system through the special point where 
the two regions of elliptic motion meet, n\ = k>2, avoids additional bifurcations. We will see 
the same behavior in the integrable billiard with parabolic boundary in a gravitational field. 



5 The gravitational billiard 

In the gravitational billiard we have a linear potential U = (g, q) with g = (0, gf. Using the 
Jacobian M of the free motion 

AJ dX(t) ( 1 tl \ 

the general solution is given by 

*«>- M (£)-(T)- (48) 

M is the same as for the ordinary billiard without potential. In this system a period one 
orbit exists if there is a point on the boundary with horizontal tangent, such that q(r) = qi 
and p(r) = — Pi- Their trace is given by 

tr[ 4) = trMK = tr(2W + rR) = r trR (49) 

and with t = (1, 0), n = (0, 1) and p = 0we obtain from (38) that 

tr^ 4) = 2t(-kti + g/n) = 4 - 2ktti. (50) 

In the latter equations we used r = 2n/g to eliminate g from the expression. With h = n 2 /2 

(2) 

and eliminating r we find tr\ = 2 — 8nh/g, which is well known, see e.g. [16, 9]. Additional 
relations like rg = 2n in this simple example, exist between the quantities entering the 
stability formulas, because we have circumvented the use of the solution of the equations 
of motion. The additional relations appear in the monodromy matrix in the S coordinate 
system: It must have the form (39), i.e. the 0's and l's are always there. In more complicated 
applications it might be more appropriate to calculate these relations directly, but the final 
check with (39) should always be made. 

The interpretation of the bifurcation of the periodic orbit is as follows: If the boundary 
is defocusing, as e.g. studied in [33], the orbit is direct hyperbolic. At a flat boundary we 
have a direct parabolic bifurcation. For a focusing boundary the orbit is elliptic for small 
energies and becomes inverse hyperbolic if the energy becomes large compared to the radius 
of curvature. The invers parabolic transition at g = 2hn typically creates an elliptic period 
two orbit in a period doubling bifurcation. Now we turn to the calculation of the stability 
of the period two orbit created in this bifurcation. We only treat the case of the time and 
space symmetric period two orbit. The spatial symmetry is a reflection at the vertical axis, 
such that the angles ccj between the tangent at the reflection points and the horizontal axis 
are related by «2 = — «i- The time reversal symmetry requires p\ = p2 = 0. 

Instead of calculating the full symmetric periodic two orbit we introduce an upright wall 
at the place of the vertical symmetry axis for which again p2 = must hold. The reflection 
matrix K2 for this wall with ot2 = ±7r/2 and K2 = is given by 
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Figure 3: Stability diagram for the symmetry reduced period two orbit in the gravitational billiard. 
The contour lines of the residue R' = 1 — R' (53) are shown. The elliptic region is shaded. Inside this 
region the lines with R' = 1/2, 3/4 are shown. In the unstable regions the integer valued residues are 
shown. The dashed line indicates the stable periodic orbit in the family of billiards in the paraboloid. 
The dotted line corresponds to the wedge billiard and the long dashed line to the circle billiard. 



Formula (42) is also valid for the gravitational billiard because it has the same Jacobian as 
the ordinary billiard, and it reduces to 

tr (4) r' t' 
2 = cos(2q ± tt) + — tr(WoRi) = — cos 2a H cos 2a(nn 2 + g cos a). (52) 

The time r' is just half the time of the full orbit. Using the relation gr' = ncosa obtained 
from the explicit solution with the above boundary values, we find 

1 - R> = R' = cos 2a{riKT' - sin 2 a) (53) 

The residue for the full orbit is now obtained from 

1-R = R=(1- 2R') 2 . (54) 

The reason that the residue (or its complement) factorize in the above fashion is usually 
related to a symmetry in the system. If we shift the zero of the potential energy to the initial 
point we can again use h = n 2 /2 to eliminate r. 

Let us now study the stability of this periodic orbit in a diagram similar to Fig. 2. Recall 
that we are not discussing the question of existence of the orbit. In fact this orbit (even in 
symmetric boundaries) does not exist for all energies. But if it exist, its stability can be read 
off Fig. 3, where the residue R' = l — R'oi the symmetry reduced orbit is plotted as a function 
of tikt' and a. We start by analyzing classes of orbits with constant angles a = 0, 7r/2,7r/4. 

For a = we have the orbit bouncing up and down on a vertical line. This is not, however, 
the period one orbit discussed above, because there is also an upper wall. The residue is just 
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1 — ukt'. For infinite energies this line describes periodic orbits of the billiard without 
potential. To compare to Fig. 2 we must pass to the unreduced system with R = 4R'(1 — R') 
and t = 2t'. Then the line a = can be mapped to the diagonal of Fig. 2. 

In the other extreme we have a = tt/2 which means a reflection at a vertical wall with 
R' = ukt'. This motion can only be realized for infinite energies or vanishing distance and 
we can connect this line to the period two orbit in the billiard without potential: the residue 
of the full system on the line a = tt/2 can again be mapped to the diagonal of Fig. 2. 

There is the surprising feature that at a = tt/4 the orbits are invers parabolic, indepen- 
dently of all the other quantities. For an intuitive explanation consider the non-symmetry 
reduced system: to throw a ball the farthest distance with least energy it must take off with 
an angle of 7r/4. For higher energies and the same distance there are two distinct solutions, 
the higher orbit taking a longer time than the lower one. Together these two orbits together 
give a period two space symmetric orbit which does not posses time reversal symmetry. In 
the full system the time symmetric orbits have residue 0, and we have a symmetry breaking 
pitch fork bifurcation, which appears as a period doubling bifurcation in the reduced system. 

Let us now look at orbits with arbitrary angle a and defocusing boundary, ukt' < 0. For 
a < 7r/4 the orbits are always inverse hyperbolic (in the reduced system). Note that in the 
full system the symmetric period two orbit is always direct hyperbolic if it is unstable. If 
the angle is increased beyond ir/4 the defocusing boundary orbit always becomes first elliptic 
and eventually direct hyperbolic as the ordinary billiard at a = tt/2 is approached. For a 
focusing boundary (ukt' > 0) the orbit is (almost) always stable at ukt' = 1/2. We observe 
that the whole picture is invariant under a rotation by tt around the point (7r/4, 1/2). With 
this duality we can translate the statement "defocusing boundary for small angles is inverse 
hyperbolic" to "sufficiently focusing boundary and large angles is inverse hyperbolic", and 
similarly "sufficiently focusing boundary and small angle is hyperbolic" . 

Now we discuss the evolution of this orbit in specific examples. The creation of the period 
two orbit by period doubling from the upright bouncing motion happens at a = 0, ukt = 1 
with residue 0. This point is part of two well known billiards, the integrable parabola billiard 
[9, 16] and the (gravitational) billiard inside the circle [15]. The line of the period two orbits 
of these systems, parametrized, e.g., by the energy, is also given in Fig. 3. In the circle billiard 
nr is constant (k is trivially a constant). For small energies (and therefore angles) the orbit 
is stable, at a = 7r/4, it becomes unstable in a period doubling bifurcation and stays unstable 
henceforth. At infinite energy the billiard becomes the integrable circle billiard without 
potential and the period two orbit is direct parabolic. For the parabola billiard the behavior 
is quite different. The creation of the orbit starts at the same place, but according to the 
already stated observation that integrable systems avoid bifurcations it stays in the elliptic 
region when the angle and energy are increased, passes to the symmetry point of the diagram 
and stops at the infinite energy limit. The third line shown in Fig. 3 corresponds to the 
period two orbit in the wedge billiard [11, 12], for which the curvature is zero. At a = tt/2 
this line joins the one of the parabola billiard. Here, however, we do not have an infinite 
energy limit but instead zero distance: The angle a is also the wedge angle, which is the only 
system parameter, and at tt/2 the wedge degenerates into a line. For a > tt/4 the orbit is 
elliptic, at 7r/4 where the billiard becomes integrable it is invers parabolic and for smaller 
angles it is inverse hyperbolic, as must be the case because the billiard is ergodic for a < tt/A 
[34]. 

Let us finally remark that the above stability formula for the full orbit is quite similar 
if the curvatures are different at the reflection points, R = (1 — 2^)(1 — 2R' 2 ), where R[ is 
given by (53) with the curvature k replaced by Kj. 
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6 Billiards with constant magnetic field or harmonic potential 



The stability formulas for the symmetric period two orbits in the billiard with harmonic 
potential or magnetic field are structurally similar to those of the ordinary billiard. In the 
symmetric case of the ordinary billiard we find from (46) that R = 2/d(l — nl/2) where I 
is the length of the orbit. The factor 1 — nl/2 will turn out to be present in the following 
examples with potential, such that its geometric significance is independent of the additional 
potential: If the Euclidean distance of reflection points is equal to the diameter of a circle 
with radius (of curvature) 1/k the stability changes. 

Constant magnetic field 

The vector potential of the constant magnetic field B is A = (—yB/2,xB/2), where we 
denoted the components of q by (x,y). The free motion takes place on circles with Larmor 
radius |v|/|B|. The Jacobian of the free motion is given by 

/ (l + D(B*))/2 -(1-D(B*))S/B 
{l> \ B(l - D(B*))S/4 (1 + D(Bt))/2 

where D(Bt) is a rotation matrix, which will always be written with its argument to dis- 
tinguish it from the rotation D around the angle a. The general solution is just given by 
X(t) = M.(t)X(0). We take B < such that the particles rotate counterclockwise. 

Since the circles of the free motion are already periodic orbits the billiard does not have 
period one orbits. For a period one orbit we need a free orbit that returns to the initial point 
in configuration space, but with different velocities. The billiad wall is then adjusted in such 
a way that the final velocity is reflected to the initial one. We again restrict ourselves to a 
symmetric period two orbit, i.e. Q2 = ol\ + 7r. Because of isotropy of space we can set ol\ = 
and «2 = vr (the same argument could have been used for the billiard without potential), such 
that W2 = Wi = Wo, where the last index denotes the standard reflection matrix (a = 
in (35)), not the position of reflection. Using the relation 2/3 = \B\t, where 2(3 is the angle of 
the part of the circle traveled by the free orbit to eliminate the time, such that p = | v| cos (3 
and n = | v| sin (3, we find 

R= (l + 2p/B/«i)(l + 2 P /Bk 2 ). (56) 

Noting that the Euclidean distance I between the reflection points is just I = \2 cos /3|v|/5| = 
\2p/B\, and that 5<0we obtain the result given in [18] for (3 < ir/2. The above formula is 
also valid for the case of the outer billiard (e.g. outside the sphere, see [19]) with (3 > ir/2. 
The stability diagram is the same as for the billiard without magnetic field given in Fig. 2. 

Harmonic oscillator 

Let the potential be given by U(q) = lj 2 x 2 /2 + u 2 y 2 /2. The Jacobian of the flow is 

M « = ( -ns, V ) • < 57 > 

with C = diag(coswt,cos^t), S; = diag(sin(jt, sin ut), and fl = diag(uj , u) . The general 
solution is X(t) = ~M.(t)X(0). For a repulsive potential are imaginary which leads to 
cosh and sinh instead of cos and sin. This system can have many different period one orbits. 
For periodicity we require q(r) = qi = (^i,yi)* / where r is the (as yet undetermined) 
period of the orbit. The condition for a period one orbit reads 

qi = C q 1 + S i r2- 1 pi- (58) 
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Excluding the special orbits with det(l — C Q ) = respectively det Si = 0, we can solve for 
qi. Inserting into the corresponding solution for p(r) gives the final momentum after time r 
again at the position qi, 

p(r) = (-SiO(l - Co)- 1 ^- 1 + C )pi = -Pi- (59) 

Because of p(r) = — pi these orbits are self retracing: they go up the potential until they 
reach the oval of zero velocity and then return the same path in configuration space until 
they reach the initial point with reversed velocity. In order to obtain a period one orbit of 
the billiard the boundary at qi must be perpendicular to pi. This is the standard situation 
for billiards without magnetic field and we already encountered it for the period one orbit in 
the linear potential. By construction these orbits have time reversal symmetry. The normal 
n at the point of reflection is parallel to pi and p = 0. Taking the trace of KM we find 

trj 4 ^ = 2 cos 2a(cos lot — cos vt) + — ^ sinioT + sin vt (60) 

LO V 

with Ru given by (36). For the calculation of R note that n is eigenvector of W and v || n, 
such that SWv = — Sv. From (58) we obtain 

LOT VT f 

Pi = (Xiuj tan — , y-yv tan — ) . (61) 

The relation between energy and period for fixed initial and final point qi is a transcendental 
equation, 

x\u? y\v 2 
2 cos 2 lot /2 ^ 2cos 2 z/r/2' { } 

and therefore we leave everything parametrized by the period r. The divergence of the energy 
as a function of r is related to the fact that the period of oscillations in the harmonic oscillator 
is independent of the amplitude. The only variation in r comes from the fact that we do 
not start at the origin. If the energy is large, the difference in r due to this initial offset is 
relatively small and for infinite energy r approaches half the period tt/lo or ir/v. 

Now we study the orbits whose reflection condition (ti,pi) = with pi given by (61) 
is fulfilled independently of r, i.e. for which x\ = or y\ = 0. With an attractive potential 
the periodic orbits with small r run upwards in the potential. For lot or vt > it the initial 
momentum in the corresponding direction changes its sign, and the orbit goes down first, 
then up on the other side and returns to qi. Periods r larger than the full period don't give 
periodic orbits in the billiard. In a repulsive potential the orbit goes up in the direction of the 
origin because otherwise it would never return. The outer billiard with repulsive potential is 
a scattering system without periodic orbits. In the case y\ = and x\ > we obtain 

(4) / LO LOT VT\ o VT i \ 

tr\ ' = 4 I 1 ± — Kx± tan — tan — J cos — , (63) 

where ± is the sign of the initial momentum. A similar expression holds for the case x\ = 0. 
In the case of the isotropic harmonic oscillator with lo = v we get the simple expression 

tr^ 4) = 4cos 2 cj^ ± 4xiKsin 2 (64) 

Introducing the energy h and U\ = U(qi) gives 

R= (l-y)(lT K 4 (65) 

In the case that h = U±, this describes a degenerate orbit which does not move at all. Let 
us first consider orbits with positive initial momentum. For larger h the orbit moves an 
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increasing distance radially outward. If the boundary is sufficiently defocusing (kx < —I) 
the orbit is always hyperbolic, while for slightly defocusing boundary with — 1 < kx < it 
is always elliptic. For a focusing boundary the orbit is stable for small energies and inverse 
hyperbolic for large energies. Note that the stability can be read off Fig. 2 with an appropriate 
reinterpretation of the axes. For orbits with initially negative x-momentum the sign of n is 
defined with respect to the opposite normal, such that both orbits with either sign of the 
initial momentum have the same stability, despite the fact that one "sees" a focusing and the 
other one a defocusing boundary. 

In a repulsive potential we have U\ < h < 0. The upper bound on h ensures that the oval 
of zero velocity around the origin (which is needed for the period one orbit) does not vanish. 
We obtain the same stability formula, however now the minus sign in front of k describes the 
(only possible) orbits with initial negative momentum. 

The symmetric period two orbits can also be easily calculated, compare to the results 
in [9]. We start at (xi,0), X\ > always with negative initial momentum until we reach 
(— x\,0). Due to the assumed symmetry we can introduce a reflecting wall with curvature 
zero in the middle and for the symmetry reduced orbit in the attractive potential we find 

# = (l-^)(l-«J/2) (66) 

and R = 4i?'(l — R 1 ) for the full orbit with I = 2x\ the distance between the reflection points. 
The factor 1 — U\/h is always positive and smaller than 1. Again the stability can be read 
off Fig. 2. For the repulsive potential we get the same formula, however, now it gives the 
complement of the residue for the symmetry reduced orbit 

# = (l-^)(l-«J/2) (67) 

and we obtain R = AR'(1 — R') for the full orbit. The existence condition in the repulsive 
case is h > 0, such that the orbit can cross the potential hill at the origin. Therefore \ — U\/h 
is always larger than 1. 

P. Stifter [35] has shown that the ordinary billiard in the cardioid, which is ergodic [5], 
is equivalent to the billiard with repulsive potential — x 2 — y 2 inside the unit circle with 
center (1, 0) for the energy H = 0. This is exactly the energy where the oval of zero velocity 
vanishes and there is an unstable equilibrium point at (0,0). Changing the energy gives a 
system which does not correspond to the cardioid billiard. Instead this gives an interesting 
example where a bounded system, which is close to integrable for low energies H ~ —4 and 
high energies, is proven to be ergodic for the special intermediate energy value H = 0. 



7 The rotating billiard 

A rotating billiard, as introduced in [20], consists of an ordinary inner or outer billiard 
boundary that moves on a circle around the origin. If we study the outer billiard with two 
discs the system shows a striking similarity to the restricted three body problem, because in 
the latter system orbits with close encounters to either large body are effectively reflected, 
even though the gravitational force is attractive, such that the hard core potential of the 
billiard can model this motion. This special case of a rotating billiard will be referred to 
as the restricted three body billiard, abbreviated as R3BB. The Hamiltonian in a uniformly 
rotating frame of reference is given by 

1 111 

H = -p 2 - u(xp y - yp x ) = -{p x + uy) 2 + -{p y - uxf - -wq 2 (68) 
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Figure 4: Stability diagram for the period one orbits in the R3BB as a function of nx and lot 
as given by the residue R in (75). The elliptic region is shaded. Inside this region the lines with 
R = 1/2, 3/4 are shown. In the unstable regions the residues —1, ±2, ±3, ±5, ±10, ±25 and ±50 are 
shown. The typical geometric situation is to the left of the vertical line kx = 1, where the boundary is 
defocusing or only slightly focusing. We consider this region now: The horizontal lines indicate saddle 
center bifurcations. The retrograde orbit created at lot w 4.08 is stable for slightly larger period, 
similarly the other retrograde orbits with lot w (2m + 1 + l/2)7r, for which h > 0. The direct orbits 
(h < 0) bifurcating at lot w (2m + l/2)ir are stable for slightly smaller period. The infinite residues at 
lot = 2m7r prevent the smooth continuation of a retrograde orbit to another retrograde orbit, similarly 
for the direct orbits. In the left two pictures the energy h and the angle 9 between in- and outgoing 
velocities at the reflection point as a function of lot are shown. 



such that we have a harmonic repulsive potential U = — cuq 2 /2 and a vector potential A = 
(—uy,ujx) t . The Jacobian of the free motion is 

M(t) _ ( DM) .DM, ) m 

such that the general solution can be written as X(t) = M(i)X(0). In the context of the 
restricted three body problem the analog of the above Hamiltonian is referred to as the 
Jacobi integral; we will stay with the former name and also call the value of the Hamiltonian 
the energy, which is not the energy of the particle in the non rotating frame. In rotating 
coordinates the discs are fixed on the x-axis at some distance from the origin. The straight 
line motion looks rather (nice and) complicated in the rotating frame of reference. The free 
motion is circle-like due to the coriolis field but it is somewhat distorted due to the repulsive 
harmonic potential centered at the origin. Depending on the value of the Hamiltonian there 
can be up to an infinite number of period one orbits coexisting. Now we are going to calculate 
the stability of these orbits. 

To have a solution that returns to qi 7^ after time r we must choose the initial momenta 
according to q q = q(r) and with (69) this gives 

qi = D(wT)qi + tD(wt)pi. (70) 

The final momenta are p(t) = D(o;r)pi and the angle between the initial and final velocities 
v = p — «4(qi) is given by 

n (1 — coswr) 2 ,„„ . 

cose = 1-2- — i J. (71) 

2 — 2 cos lot — 2lot sin lot + lo^t^ 
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For lot = 27rm, m = ±1, ±2, . . ., is zero such that the free motion is periodic by itself; 
in the non rotating frame the particle does not move at all. This is similar to the ordinary 
billiard, where an orbit touching the boundary becomes infinitely unstable due to n — > 0. 
For every other value of r, however, a periodic orbit of the billiard exists. The angle between 
the normal at the reflection point and the velocity in the non rotating frame (!) is given by 
i9 = (ir — lot)/2 mod it — tt/2. For negative # the velocity is directed against the direction 
of rotation of the disk (retrograde), while for positive # the direction is the same (direct). 
Retrograde orbits typically have higher energy than direct orbits because they have to be 
faster to reach the next collision. The energy for these periodic orbits parametrized by r is 
obtained by eliminating p in (68) using (70), 

h = —(I — cos lot — lot sin lot), (72) 

which is shown in Fig. 4, where |qi| = r = 2. Every intersection of the line h = const with 
this curve corresponds to a periodic orbit of the system. At extrema of this curve periodic 
orbits are created/destroyed in pairs when the energy is varied. At maxima retrograde orbits 
bifurcate and at the minima direct orbits bifurcate. These important points are given by 

h' = — = o (2 — 2 cos lot — 2lot sin lot + lo 2 t 2 cos lot) = 0, (73) 

OT T 6 

which will turn out to be closely related to the stability. 

Taking the trace of KM we find the general stability formula to be 

tr[ 4) = r {coslot{Ru + R 22 ) + sinuoT{R 2 i - Rut)) , (74) 

and for the use of (38) we need (d„U,n) = L0 2 (xisina — yicosa) and B = —2lo. In order 
for the reflection to map v(r) into v(0) we must have (p(r) — „4(qi),ti) = ((pi — „4(qi),ti), 
which gives a.\ = arctan(x 1 /y 1 ). We now choose a coordinate system in such a way that the 
orbit starts at y\ = 0, such that a\ = tt/2 and D = S. Introducing the distance from the 
origin r = x\ the residue obtained from (74) is 

12 — 2cos lot — 2lot sin LOT + L0 2 T 2 COS LOT . . , 

R = (r/e-1) (75) 

2 cosu;t-1 v ' y ' 

, 9 9 COS LOT — l w 

= (^^-^-pX^-l) (76) 

= . (77 ) 

2r z 1 — cos lot 

The last formula nicely illustrates the saddle-center bifurcations at the extrema of H(t): Two 
new orbits are created with stability R = 0, i.e. h! = 0, which is just the condition for the 
bifurcation. Close to the bifurcation point h! has different signs for the two orbits, such 
that one is elliptic and the other one direct hyperbolic. The elliptic orbit becomes inverse 
hyperbolic typically after only a slight change of r. The interval of stability decreases with 
increasing r. The whole scenario is illustrated in the stability diagram of Fig. 4, which 
completes the results reported in [20]. 

As our final application we now turn to the calculation of the period two orbits in the 
R3BB, which are of great interest because of their apparent similarity with orbits in the 
restricted three body problem that have close encounters with both centers. We treat the 
symmetric case with tt = Ki and qi = (r, 0) and q2 = (— r, 0). The analogue of (70) gives pi 
for the given qi and q2- The resulting energy is 

- r 2 

h = —(I + cos lot + lot sin lot). (78) 
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The main difference to the period one orbits is that h diverges for r — ► 0. Therefore we have a 
period two orbit up to infinite energies, which bounces back and forth between the two disks 
with the motion becoming a straight line for r — > and h — > oo. Except for this behavior we 
again have the oscillations in h(r) and therefore a bifurcation scenario similar to that of the 
period one orbits. The important features of the residue again depend on the derivative of 
h, which is given by 

dh r 2 

h! = — — = o (2 + 2cosc<jr + 2ujt smujT — uj 2 t 2 cos ujt). (79) 

Because of the above symmetry assumptions the t-n-coordinates are given by ol\ = tt/2 and 
«2 = 37r/2 such that Ki = K2 with Wi = W2 = Wo- Although the orbit is not time reversal 
symmetric we have r = t\ = T2 due to the spatial symmetry and therefore Mi = M2. For 
the trace of KMKM we find 

t4 4) =4(2#-l) 2 , R' = ~ti^ ™ - 1 , (80) 
2 v ' 2r 2 1 + coswr v ' 

giving the same factorization of R = ti^f' '/4 as in the gravitational billiard (54). 

The bifurcation scenario for the period two orbits is a follows. They are created in saddle 
center bifurcations at the extrema of h. The residue of the elliptic orbit increases from 
to 1. However the residue does not increase beyond 1 but instead begins to decrease again. 
This "touching of residue 1" is an effect of the symmetry of the orbit, which corresponds to 
passage of R' through 1/2. There are two period four orbits involved in this 2/4-bifurcation 
[30] and the old orbit stays elliptic. Eventually it reaches R = again and looses stability 
in a period doubling bifurcation. For retrograde orbits this scenario happens for decreasing 
period r, and for direct orbits with increasing period. Decreasing the period of the retrograde 
orbits further eventually turns them infinitely unstable at 1 + cos ujt = 0, where they have a 
tangency instead of a reflection. The same happens for the direct orbits if we increase their 
period. However, if we increase the period of a retrograde orbit it becomes more and more 
unstable, but eventually becomes direct and then stable again, thus every two of the orbits 
described above are continuously connected in phase space. The whole scenario is repeated 
for any number m of revolutions that are completed in the non rotating frame before the 
next reflection. The intervals of stability are extremely small for higher periods. Only in 
the case nr « 1, i.e. the distance from the origin is close to the radius of curvature, are the 
windows of stability larger. For kt = 1 all the orbits become inverse parabolic, independently 
of the period. Note that this case is completely different from the standard case in the R3BB, 
because the curvature must be positive. As a matter of fact, the case kt = 1 can be realized 
as an orbit inside a circle, and we reobtain the integrable circle billiard described in rotating 
coordinates. The distance between the reflection points of the orbit is I = 2r, such that the 
same geometric factor enters the residue as in the previous examples. 

It is possible to treat the general period two orbit, with different curvatures and distances 
from the origin. This case is quite promising in order to make a more quantitative comparison 
with the restricted three body problem, because with these parameters we have an analog 
of the mass ratio in our system. In any case the system at hand, which was introduced in 
[20], gives the opportunity to obtain analytical stability results. This tempers the feeling 
that this approximation of the three body problem is too crude, and we hope, as cited in the 
beginning, that "only the interesting qualitative questions need to be considered" . 

8 Conclusion 

We have derived a general expression for the linearized Poincare map in a billiard with 
potential. Due to the additional coupling introduced by a (vector) potential the linearized 
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map can only conveniently be written in four dimensions. Using a canonical transformation in 
phase space we have calculated the linearization and it turned out that the arc length and the 
parallel component of the velocity (not of the momentum) give a canonical coordinate system 
on the surface of section given by the billiard boundary. This holds true for any potential. 
The resulting formula is useful for the application of numerical methods to find periodic 
orbits, or for the estimation of the Lyapunov exponent. Moreover, since the linearization in 
phase space includes the variation of the energy, which usually gives a nontrivial change in 
the orbit for a billiard with potential, it can also be used to follow periodic orbits when the 
energy is varied. 

The main application of the linearized map is the calculation of the stability of periodic 
orbits. The trace of the (four dimensional) monodromy matrix of the periodic orbit can be 
factorized into a product of matrices describing the piecewise free motion between reflections 
and the contribution from the reflection. This generalizes a well known formula with 2x2 
matrices for ordinary billiards. With these results it was possible to perform the stability 
analysis of period one and two orbits in four billiards with different (vector) potentials. 

The results are presented in form of stability diagrams, where the residue of the orbits is 
shown in its essential parameter dependence. These diagrams can give some intuition about 
the mechanisms that create stability respectively instability. The investigation of stability 
diagrams of the ordinary billiard (which is the same as for the one with magnetic field) and 
the gravitational billiard showed that periodic orbits in one parameter families of integrable 
systems tend to avoid period doubling bifurcations which are generically present. In the 
final application to a rotating billiard, which mimics the restricted three body problem, we 
have obtained an interesting bifurcation scenario of an infinite number of period one and two 
orbits by successive saddle-center bifurcations with increasing period. The calculations for 
this system can be extended to the non-symmetric case, in order to obtain a simple system 
that presents some essential features of the restricted three body problem in an analytically 
tractable way. 
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